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Fourier normal ordering 36 is a new algorithm to construct explicit rough 
paths over arbitrary Holder-continuous multidimensional paths. We apply in 
this article the Fourier normal ordering ordering algorithm to the construc- 
tion of an explicit rough path over multi-dimensional fractional Brownian 
motion B with arbitrary Hurst index a (in particular, for a < 1/4, which 
was till now an open problem) by regularizing the iterated integrals of the 
analytic approximation of B defined in [34j . The regularization procedure is 
applied to 'Fourier normal ordered' iterated integrals obtained by permut- 
ing the order of integration so that iimermost integrals have highest Fourier 
modes. The algebraic properties of this rough path are best understood us- 
ing two Hopf algebras; the Hopf algebra of decorated rooted trees jB] for the 
multiplicative or Chen property, and the shuffle algebra for the geometric or 
shufHe property. The rough path lives in Gaussian chaos of integer orders 
and is shown to have finite moments. 

As well-known, the construction of a rough path is the key to defining a 
stochastic calculus and solve stochastic differential equations driven by B. 

The article [37] gives a quick overview of the method. 
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Introduction 



The (two-sided) fractional Brownian motion t ^ Bt, t G M (fBm for short) 
with Hurst exponent a, a G (0, 1), defined as the centered Gaussian process 
with covariance 

E[BsBt] = ^(|sp" + -\t- s|2"), (0.1) 

is a natural generalization in the class of Gaussian processes of the usual 
Brownian motion (which is the case a = ^), in the sense that it exhibits 
two fundamental properties shared with Brownian motion, namely, it has 
stationary increments, viz. E[{Bt-Bs){Bu-By)] =E[{Bt+a-Bs+a)iBu+a- 
By^a)] for every a,s,t,u,v E M, and it is self-similar, viz. 

VA > 0, {Bxt,t G M) ^ (X^Bt, t E M). (0.2) 

One may also define a d-dimensional vector Gaussian process (called: d- 
dimensional fractional Brownian motion) by setting Bt = {Bt{l), . . . , Bt{d)) 
where {Bt{i),t G IK)i=i,...,d are d independent (scalar) fractional Brownian 
motions. 

Its theoretical interest lies in particular in the fact that it is (up to 
normalization) the only Gaussian process satisfying these two properties. 

A standard application of Kolmogorov's theorem shows that fBm has 
a version with a~ -Holder continuous (i.e. k- Holder continuous for every 
K < a) paths. In particular, fBm with small Hurst parameter a is a natural, 
simple model for continuous but very irregular processes. 

There has been a widespread interest during the past ten years in con- 
structing a stochastic integration theory with respect to fBm and solv- 
ing stochastic differential equations driven by fBm, see for instance \24:\ 
\TE\ [5l [30l I31j . The multi-dimensional case is very different from the one- 
dimensional case. When one tries to integrate for instance a stochastic 
differential equation driven by a two-dimensional fBm B = (B{1), B{2)) by 
using any kind of Picard iteration scheme, one encounters very soon the 
problem of defining the Levy area of B which is the antisymmetric part of 
Ats := £dBt,{l) dBt^{2). This is the simplest occurrence of iterated 

integrals B^^{ii, . . . ,ik) := Jl dBt^{ii) . . . Jl'"'^ dBt^{ik) , ii,...,ik < d for 
d-dimensional fBm B = {B{1), . . . , B{d)) which lie at the heart of the rough 
path theory due to T. Lyons, see [25^ I26j. An alternative construction has 
been given by M. Gubinelli in [T6] under the name of 'algebraic rough path 
theory', which we now propose to describe briefly. 
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Assume Tt = (rt(l), . . . ,Tt{d)) is some non-smooth d-dimensional path 
which is a-Holder continuous. Integrals such as j fi{Tt)drt{l) + ... + 
fd{Tt)dTt{d) do not make sense a priori because T is not difFerentiable 
(Young's integral [23j works for a > | but not beyond). In order to de- 
fine the integration of a differential form along F, it is enough to define a 
geometric rough path (T^ , . . . ,T^^^"^) lying above F, [l/aj=entire part of 
1/a, where = (5F)ts := F^ — Fg is the increment of F between s and 
t, and each F'^ = (F'^(«i, . . . , ifc))i<j^^, ^jj.<(^, /c > 2 is a substitute for the 
iterated integrals /J dFt^ (ii) J^^ dTt2 (^2) • • • Js''~^ dVt^ {ik) with the following 
three properties: 

(i) {Holder continuity) each component ofT^ is ka~ -Holder continuous, 
that is to say, kn-Hdlder for every k < a; 

(a) (multiplicativity) letting 5V\^g := T^^ — T^^ — F^^, one requires 

5FL(n,...,4)= r,^^(n,...,4jr^K4,+i,...,4); (o.s) 

ki+k2=k 

(Hi) (geometricity) 

^4/(^1, • • • j^ni)rjs^(jl, . . . , jn2) = r^/"'' ^(fci, . . . , /c^^+nj) 

fc6Sh(i,j) 

(0.4) 

where Sh(i,_7) is the subset of permutations of ii, . . . ,ini, ji, ■ ■ ■ , jn2 
which do not change the orderings 0/ (ii, . . . , i„^) and (ji, . . . , Jn2)- 

The multiplicativity property implies in particular the following identity 
for the (non anti-symmetrized) Levy area: 

Ats = Atu + Aus + {Bt{l) - B.,{l)){Bu{2) - 5,(2)) (0.5) 

while the geometric property implies 

/ dBt,{l) / dBt2{2)+ / dBt2{2) / dBt,{l) 

•J s ^ s ^ s ^ s 

fdBt,{l)^ (£di?i,(2)) = {Bt{l) - Bs{mBt{2) - B,{2)). 

(0.6) 

Then there is a standard procedure which allows to define out of these 
data iterated integrals of any order and to solve differential equations driven 
by F. 



4 



The multiplicativity property (j0.3p and the geometric property (|0.4p are 
satisfied by smooth paths, as can be checked by direct computation. So the 
most natural way to construct such a multiphcative functional is to start 
from some smooth approximation T^, rj ^ of T such that each iterated 
integral r^J^(zi, . . . , i^), k < [1/aJ converges in the A;K-Holder norm for 
every k < a. 

This general scheme has been applied to fBm in a paper by L. Coutin and 
Z. Qian [9j and later in a paper by the author [33], using different schemes 
of approximation of B by with r] ^ 0. In both cases, the variance of the 
Levy area has been proved to diverge in the limit r/ — > when a < 1/4. 

The approach developed in [34j makes use of a complex-analytic process 
r defined on the upper half-plane 11+ = {z = x + iy \ y > 0}, called F- 
process or better analytic fractional Brownian motion (afBm for short) |33j . 
Fractional Brownian motion Bt appears as the real part of the boundary 
value of when Im z ^ 0. A natural approximation of Bt is then obtained 
by considering 

Bl := Tt+ir^ + Tvi^ = 2Re Tf+i^^ (0.7) 

for r] ^ 0. We show in subsection 3.1 that B^ may be written as a Fourier 
integral, 

B^ = cJ e-^l«l|e|^-°^^^ W{dO (0.8) 

for some constant Cq,, where {W{^),^ > 0) is a standard complex Brownian 
motion extended to M by setting W{—^) = —W{^), ^ > 0. When ^ 0, 
one retrieves the well-known harmonizable representation of B |32j . 
The so-called analytic iterated integrals 

' h{zi)dTM r h{z2)dT,,{2) . . . fd{zd)dT,^{d) 

J s J s 

(where fi, ■ ■ ■ , fd are analytic functions), defined a priori for s,t £ 11+ by 
integrating over complex paths wholly contained in 11+ , converge to a finite 
limit when Im s,Im t — > |34j . which is the starting point for the construc- 
tion of a rough path associated to F [33]. The main tool for proving this 
kind of results is analytic continuation. 

Computing iterated integrals associated to = 2 lim^^o Re F^+i^ in- 
stead of F yields analytic iterated integrals, together with mixed integrals 
such as for instance dF^-^ (1) f^^ (iF22(2). For these the analytic continu- 
ation method may no longer be applied because Cauchy's formula fails to 
hold, and the above quantities may be shown to diverge when Re s. Re t ^ 
0, see [SI [35]. 
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Let us explain first how to define a Levy area for B. Proofs (as well 
as a sketch of the Fourier normal ordering method for general iterated in- 
tegrals) may be found in [37]. As mentioned before, the uncorrected area 
Al := Jl dBl^{l) J^^ dB'l^{2) diverges when r? ^ 0+. The idea is now to 
find some increment counterterm {6Z'^)ts = — Zg such that the regu- 
larized area IZA^g ■= A^^ — {5Z'^)ts converges when t/ ^ 0+. Note that 
the multiplicativity property (jO.Sp holds for TZA^ as well as for A^ since 
{5Z^)ts = {6Z^)tu + {SZ^)us- This counterterm Z^ may be found by using 
a suitable decomposition of A^g into the sum of : 

- an increment term, {5G^)ts', 

- a boundary term denoted by A^g{d). 

The simplest idea one could think of would be to set 

i5G^)ts= fdBl{l)Bli2), (0.9) 

J s 

and 

Alid) = - fdBlil) . B^{2) = -B^,{2){Bni) - B^,{1)). (0.10) 

J s 

Alternatively, rewriting A^g as dBu2{2) J^^ dBui{l), one may equivalently 
set 

{SG%s = - fdBl{2)Bl{l) (0.11) 

J s 

and 

Alid) = fdBl{2) . ^,"(1) = i?I'(l)(i??(2) - 5,^(2)). (0.12) 

J s 

Now (5G^ diverges when r/ ^ 0^, but since it is an increment, it may be 
discarded (i.e. it might be used as a counterterm). The problem is, AKd) 
converges when — > 0^ in the K-Holder norm for every k < a, but not in 
the 2ft;-Holder norm (which is of course well-known and may be seen as the 
starting point for rough path theory). 

It turns out that a slight adaptation of this poor idea gives the solution. 
Decompose A^ into a double integral in the Fourier coordinates ^1,^2 using 
(jO.Sp . Use the first increment /boundary decomposition (j0.9|0.10p for all 
indices |6I < 161, and the second one (I0.11l0.12p if |Ci| > 161- Then AUd), 
defined as the sum of two contributions, one coming from (jO.lOp and the 
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other from ()0.12p . does converge in the 2K-Holder norm when t/ — > 0"*", for 
every k < a. 

As for the increment term 6G^ , defined similarly as the sum of two con- 
tributions coming from (j0.9p and (jO.lip . it diverges as soon as a < 1/4, but 
may be discarded at will. Actually we use in this article a minimal regu- 
larization scheme: only the close-to-diagonal (i.e. ^1/^2 ~ —1) terms in the 
double integral defining 6G^ make it diverge. Summing over an appropriate 
subset, e.g. —.^1 [^2/2,2^2] yields an increment which converges (for every 
a G (0, ^)) when — > in the 2/t-Holder norm for every k < a. 

Let a < 1/4. As noted in |35j . the uncorrected Levy area of the 
regularized process converges in law to a Brownian motion when 77 — > 0^ 
after a rescaling by the factor 7/2*-^"^"^. In the latter article, the following 
question was raised: is it possible to define a counterterm living on the 
same probability space as fBm, such that (i) the rescaled process r/2 (^"^"^X^ 
converges in law to Brownian motion; (ii) {B^,A^ — X^) is a multiplicative 
or almost multiplicative functional in the sense of [23], Definition 7.1; (iii) 

— X^ converges in the 2K-Holder norm for every k < a when r/ — > ? 
The counterterm X"^ := A"^ — TZA^ gives a solution to this problem. 

The above ideas have a suitable generalization to iterated integrals 
J dB{ii) . . . J dB{in) of order n > 3. There is one more difficulty though: 

decomposing {B'^)'^^{ij) into / dWg^. (ij)e^"J^Je~''l^-'l|^j|^~", an extension 
of the first increment/boundary decomposition (|0.91 10.10|) . together with a 
suitable regularization scheme, yield the correct Holder estimate provided 
ICil < • • • < \in\- What should one do then if < . . . < |Co-(n)l foi' some 

permutation a instead ? The idea is to permute the order of integration by 
using Fubini's theorem, and write dBu^ (ii) . . . dBu„{in) as some it- 

erated tree integral J diJ^i (V(i)) • • • / dBu„{i„(^n-^). The integration domain, 
in the general case, becomes a little involved, and necessitates the introduc- 
tion of combinatorial tools on trees, such as admissible cuts for instance. 
The underlying structures are those of the Hopf algebra of decorated rooted 
trees [61 [7] (as already noted in [21] or [T^ ) , and of the Hopf shuffle algebra 
[27l [28] . The proof of the multiplicative and of the geometric properties 
for the regularized rough path, as well as the Hopf algebraic reinterpreta- 
tion, are to be found in [36]. The general idea (see subsection 2.5 for more 
details) is that the fundamental objects are skeleton integrals (a particular 
type of tree integrals) defined in subsection 2.1, and that any regulariza- 
tion of the skeleton integrals (possibly even trivial) yielding finite quantities 
with the correct Holder regularity produces a regularized rough path, which 
implies a large degree of arbitrariness in the definition. The idea of can- 
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celling singularities by building iteratively counterterms, originated from 
the Bogolioubov-Hepp-Parasiuk-Zimmermann (BPHZ) procedure for renor- 
malizing Feynmann diagrams in quantum field theory [18j, mathematically 
formalized in terms of Hopf algebras by A. Connes and D. Kreimer, has been 
applied during the last decade in a variety of contexts ranging from numer- 
ical methods to quantum chromodynamics or multi-zeta functions, see for 
instance \21\ [28l [38] . We plan to such a (less arbitrary) construction in the 
near future (see discussion at the end of subsection 2.5). 



The main result of the paper may be stated as follows. 

Theorem 0.1 Let B = (-B(l), . . . , B{d)) he a d- dimensional fBm of Hurst 
index a G (0, 1), defined via the harmonizable representation, with the associ- 
ated family of approximations B^ , t] > living in the same probability space, 
see eq. (E^. Then there exists a rough path (7^B^''' = SB"^, . . . ,7^BLV"J.^) 
over B^ [r] > 0), living in the chaos of order 1, . . . , [1/aJ of B, satisfying 
properties (ii) (multiplicative property) and (Hi) (geometric property) of the 
Introduction, together with the following estimates: 

(uniform Holder estimate) There exists a constant C > such that, for every s,t G M and r] > 0, 

E|7^B^;''(n,...,i„)|2<C7|^-s|2""; 

(rate of convergence) there exists a constant C > such that, for every s,t gM. and r]i,rj2 > 
0, 

E|7^B^;''^(n, ...,in)- TZB^f'iii, . . .,in)\' < C\vi - r?2p". 

These results imply the existence of an explicit rough path TZB over B, 
obtained as the limit of TZB'^ when r/ — > 0. 

Here is an outline of the article. We first recall briefly some definitions 
and preliminary results on algebraic rough path theory in Section 1, which 
show in particular that Theorem 10.11 implies the convergence of TZB^ to a 
rough path T^B over fractional Brownian motion B when — > 0. Section 2 
is dedicated to tree combinatorics and to the introduction of quite general 
regularization schemes for the iterated integrals of an arbitrary smooth path 
r. The proof of the multiplicative and geometric properties are to be found 
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in [36] and are not reproduced here. We apply a suitable regularization 
scheme to the construction of the regularized rough path TZB^ in section 3, 
and prove the Holder and rate of convergence estimates of Theorem 10.11 for 
the iterated integrals 7^B"'''^(zi, . . . with distinct indices, i\ ^ . . . ^ in- 
We conclude in Section 4 by showing how to extend these results to coincid- 
ing indices, and introducing a new, real-valued, two-dimensional Gaussian 
process which we call two-dimensional antisymmetric fractional Brownian 
motion, to which the above construction extends naturally. 

Notations. The group of permutations of {1, . . . ,n} will be denoted 
by S„. The Fourier transform is : / — > Tf{i) = f f{x)e~^^^dx. If 
\a\ < C\b\ for some constant C (a and b depending on some arbitrary set of 
parameters), then we shall write \a\ < \b\. 

1 The analysis of rough paths 

The present section will be very sketchy since the objects and results needed 
in this work have ahead been presented in great details in [33]. The fun- 
dational paper on the subject of algebraic rough path theory is due to M. 
Gubinelli [M], see also [Ej for more details in the case a < 1/3. Let us re- 
call briefly the original problem motivating the introduction of rough paths. 
Let r : R ^ be some fixed irregular (i.e. not differentiable) path, say 
K-Holder, and / : M ^ R'^ some function which is also irregular (mainly 
because one wants to consider functions / obtained as a composition g oT 
where g : M."^ ^ is regular). Can one define the integral J fxdVx ? The 
answer depends on the Holder regularity of / and T. Assuming / is 7- 
Holder with k + j > 1, then one may define the so-called Young integral [23] 
Jl fxdT^ as the Riemann sum type limit \vco.\n\^oY.{tj}(in fui^U+i - ^u), 
where Yl = {s = t^ < . . . < tn = t} \s a, partition of [s,t] with mesh |n| 
going to 0. Then the resulting path Yt — Ys := fxd^x has the same reg- 
ularity as r. If K -|- 7 < 1 instead, this is no more possible in general. One 
way out of this problem, giving at the same time a coherent way to solve 
differential equations driven by F, is to define a class of T- controlled paths 
Q, such that the above integration problem may be solved uniquely in this 
class by a formula generalizing the above Riemann sums, in which formal 
iterated integrals r"'(ii, . . . , z„) of F appear as in the Introduction. 

Definition 1.1 (Holder spaces) Let k S (0, 1) and T > fixed. 

(i) Let Cf = Cf([0, r],C) be the space of complex-valued n- Holder func- 
tions f in one variable with (semi-)norm ||/||k = supg^^g[o,T] ^'^^]^r^p^^ ■ 
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(ii) Let 6*2 = C2([0,T],C) be the space of complex-valued functions f = 
fti,t2 ^"^^ variables vanishing on the diagonal ti = t2, such that 
|2,K < c>o, where \ \ . ||2,k is the following norm: 

2a = sup -. (1.1) 

s,te[o,T] |t - sr 

(Hi) Let C3 = C3([0,T],C) be the space of complex-valued functions f = 
fti,t2,t3 of three variables vanishing on the subset {ti = t2} U {t2 = 
ts} U {ti = t^}, such that ||/||3,k < 00 for some generalized Holder 
semi-norm \ \ . Ha^^ defined for instance in fl6f . section 2.1. 



Definition 1.2 (increments) (i) Let f be a function of one variable: 
then the increment of f, denoted by df, is {5f)ts ■= fit) — f{s). 

(ii) Let f = fts be a function of two variables: then we define 

{Sf)tus ■= fts — ftu — fus- (1-2) 
Note that 6 o S{f) = if f is a function of one variable. 

Let r = (r(l), . . . , T{d)) : [0, T] ^ M"^ be a /t-Holder path, and (Tliii) := 

^tih) - rs(n),rL(^i;^2),---,r[y''J(ii,...,iLi/Kj)) be a rough path lying 
above F, satisfying properties (i) (Holder property), (ii) (multiplicativity 
property) and (iii) (geometricity property) of the Introduction. 

Definition 1.3 (controlled paths) Let z = {z{l), . . . , z{d)) G Cf for some 
K < a and N = [1/kJ +1. Then z is called a (T-) controlled path if its in- 
crements can be decomposed into 

N 

'^^W=E E r"(ii,...,ij.r(ii,...,i„;i)+<70(i) (1.3) 

for some remainders g^{i) G C'^'' and some paths /"(zi,... ,in',i) G (Cf)" 
such that 

(5/"(zi,...,z„;z) = 

N-l-n 

E E ■ ■ ■ .jl)-f^'^{h, ■ ■ ■ ,jhh, ■ ■ ■ ,in]i) + g'^iii, ■ ■ ■ ,in]i] 
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for some remainder terms g^{ii, . . . ,in',i) G 

We denote by the space of all such paths, and by Q^- the intersection 

We may now state the main result. 

Proposition 1.4 (see [17j, Theorem 8.5, or [33], Proposition 3.1) Let 

z G Qcf- . Then the limit 

n d 



nt n a 

/ ZxdVx := lim 



fc=0 i=l 



N-1 

^xt,+^,t,{i)zt,{i) + Y, J2 rr,tlA(^i>--->^n,i)cr,(n,...,in;i) 

n=l 

(1.5) 

existe in the space Q^- ■ 

Assume F is a centered Gaussian process, and a family of Gaus- 
sian approximations of T living in its first chaos. Then the Proposition 
below gives very convenient moment conditions for a family of rough paths 
(r'',r^'^, . . . , pL^/'^-l'^) to converge in the right Holder norms when 77 — > 0, 
thereby defining a rough path above T. 

Proposition 1.5 Let T be a d-dimensional centered Gaussian process ad- 
mitting a version with a.s. a~ -Holder paths. Let N = [l/a\. Assume: 

1. there exists a family , 77 — > 0^ of Gaussian processes living in the 
first chaos of T and an overall constant C such that 

(i) 

Elr^^-r^P <C|t-s|2"; (1.6) 

(ii) 

E|r?-rf|2<C|e-ry|2"; (1.7) 

(Hi) yt G [0,T], Fj ^ Tt when r] 0; 

2. there exists a truncated multiplicative functional (vlf = T^—T^, T^f, . . . , T^^^) 
lying above and living in the n-th chaos of T, n = 1 . . . ,N, such 

that, for every 2 < k < N , 

(^) 
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(ii) 

E|r,^/-r,^f|2<C|e-^r. (1.9) 

Then (r^'^, . . . ,r^'^) converges inL'^{Vl- C^([0, T], M'^) xC72'^([0, T], M'^') x 
. . . X C^'^([0,r],M'^'^)) for every k < a to a rough path (r\ . . . ,r^) lying 
above T. 

Short proof (see |33j . Lemma 5.1, Lemma 5.2 and Prop. 5.4). The 
main ingredient is the Garsia-Rodemich-Rumsey (GRR for short) lemma 
[B] which states that, if / G C^{[0,T],C), 

I I/lb,. < C (^1 15/113,. + (/"^/^ iJ-Zfl^^ ^")'^'') ^'-''^ 
for every p > 1. 

Then properties (|1.6ll.8p imply by using the GRR lemma for p large 
enough, Jensen's inequality and the equivalence of L^-norms for processes 
living in a fixed Gaussian chaos 

E||r'='''||2,fc.<E||<5r^'^||3,fc. + c. (1.11) 

By using the multiplicative property (ii) in the Introduction and induction 
on k, 'E\\ST^''^\\s^kK may in the same way be proved to be bounded by a 
constant. 

On the other hand, properties ()1.6|1.7|1.8|1.9I) . together with the equiv- 
alence of L^-norms, imply (for every k < a) 

E\T'lf - r,^;^|2 <\t- s\^''^\e - r/|2("-'^) (1.12) 

hence, by the same arguments, 

E||r'='^-r^'^||2,fc.<k-^?r-'^ (1.13) 

which shows that F is a Cauchy sequence in 

2 Tree combinatorics and the Fourier normal or- 
dering method 

2.1 Prom iterated integrals to trees 

It was noted already long time ago [3j that iterated integrals could be en- 
coded by trees. This remark has been exploited in connection with the 
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construction of the rough path solution of (partial, stochastic) differential 
equations in [T7]. The correspondence between trees and itegrated integrals 
goes simply as follows. 

Definition 2.1 A decorated rooted tree (to be drawn growing up) is a finite 
tree with a distinguished vertex called root and edges oriented downwards 
(i.e. directed towards the root), such that every vertex wears an integer 
label. 

If T is a decorated rooted tree, we let V{Y) be the set of its vertices 
(including the root), and I : y(T) — > N 6e its vertex labeling. 

More generally, a decorated rooted forest is a finite set of decorated rooted 
trees. If T = {Ti, . . . ,Ti} is a forest, then we shall write T as the formal 
commutative product Ti . . . T/. 

Definition 2.2 Let T be a decorated rooted tree. 

• Letting v,w £ V{T), we say that v connects directly to w, and write 
V ^ w or equivalently w = , if {v,w) is an edge oriented downwards 
from V to w. (Note that v~ exists and is unique except ifv is the root). 

• If Vm — > Vm-i ■ . . ^ vi, then we shall write Vm vi, and say that 
Vm connects to vi. By definition, all vertices (except the root) connect 
to the root. 

• Let {vi, . . . ,f|v(T)|) be an ordering ofV{T). Assume that {vi Vj) ^ 
{i > j) (in particular, vi is the root). Then we shall say that the 
ordering is compatible with the tree partial ordering defined by 

Definition 2.3 (i) Let T = (T{1), . . . ,T{d)) be a d-dimensional smooth 
path, and T a decorated rooted tree such that i : V{T) {1, . . . ,d}. 
Then /T(r) : ^ M is the iterated integral defined as 



partial ordering. 

In particular, if T is a trunk tree with n vertices (see Fig. CP - so that 
the tree ordering is total - we shall write 



where {vi 



[hints 




(^(^^|y(T)|)) 



/T(r) = liiT) 



(2.2) 
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where 

dFM^)) / dr,,m))... / dF^Mn)). (2.3) 

(ii) (generalization) Assume T is a subtree ofT. Let n he a Borel measure 
onM^. Then 

[h^ts := f ... r>^'>l (2.4) 

J s J s J s 

is a measure on M'^^'^. 

" T 

:1 

Figure 1: Trunk tree. 

Assume T = T so [If{ij)\ts is a number. Then case (i) may be seen 
as a particular case of case (ii) with ^ = dl'{i{vi)) (55 ... '55 dV {i{v\y . 
Conversely, case (ii) may be seen as a multilinear extension of case (i), and 
will turn out to be useful later on for the regularization procedure. Note 
however that (i) uses the labels of T while (ii) doesn't. 

The above correspondence extends by (multi)linearity to the algebra of 
decorated rooted trees which we shall now introduce. 

Definition 2.4 (algebra of decorated rooted trees) (i) LetT he the 
free commutative algebra over M generated by decorated rooted trees. 
// Ti, T2, . . . T; are (decorated rooted) trees, then the product Ti . . . T; 
is the forest with connected components Ti, . . . , T;. 

(ii) Let T' = Yli=i ''^l^l ^ > where G M and each = T;^i . . . 

is a forest with labels in the set {l,...,d}, and F be a smooth d- 
dimensional path as above. Then 

L 

[iTir)]ts ■.= Y,rni[iT,Ants . . . [htMi)(r)]ts. (2.5) 
1=1 
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Let us now rewrite these iterated integrals by using Fourier transform. 

Definition 2.5 (formal integral) Let / : M — > M 6e a smooth, compactly 
supported function such that J-f{0) = 0. Then the formal integral f^f = 
-JJoffzs defined as ^ f^^ {J^f m^^ d^. 

Formally one may write: 

ft pt iti 

/ e'''^dx= / e'^'«(ix = V (2.6) 
J J±ioo i? 

(depending on the sign of ^). The condition J^f{0) = prevents possible 
infra-red divergence when — > 0. 

The skeleton integrals defined below must be understood in a formal 
sense because of the possible infra-red divergences. 

Definition 2.6 (skeleton integrals) (i) Let T be a tree with £ : T — > 
{1, . . . ,d} and T be a d- dimensional compactly supported, smooth path. 
Let {vi, . . . jV^y^j-j^) be any ordering ofV{T) compatible with the tree 
partial ordering. Then the skeleton integral of T along T is by defini- 
tion 

[Skh{r)]t = f dr,^^{i{v,)) dr,Mv2))...p™ dv^^^^^^^^iiiv^y^^^^)). 

(2.7) 

(a) (multilinear extension, see Definition \2.3\) Assume T is a subtree ofT, 
and fi a compactly supported Borel measure on . Then 

[Sk/Tr(/i)]t = f p' . . . I""'"*^" ^l{dx,, , . . . , ) (2.8) 

is a measure on M''^'^^. 



Formally again, [Sk/']r(r)]i may be seen as [/T(r)]f,±ioo- Note that (de- 
noting by fi the partial Fourier transform of /x with respect to {xv)v<^v(T))i 
the following equation holds, 



[SklT(;u)]t = (27r)- 



-|V(T)|/2. 



(/i, SkiT ( {xy 



ev(T) 



(2.9) 
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Lemma 2.7 The following formula holds: 

(2.10) 

Proof. We use induction on |T^(T)|. After stripping the root of T 
(denoted by 0) there remains a forest T' = T'^^ . . . T'j, whose roots are the 
vertices directly connected to 0. Assume 

[Sk/T;.(r)k =/■••/ n d^v.e^°^^'^^^'^^'''FMv).er,). (2-11) 



Note that 



( nsk/T;.(r) ) (0 = / _ n d^vll^Mv)vmr,))- 

\j=l J •^^vevi'i)\{0}iv-i veV{T)\{Q} j=l 



Then 



(2.12) 

[Sk/T(r)]t = r drM0))tl[Skir.{r)U 



Mia 

-oo 



(£(o))nsk/T^(r)l (0 

1 r+oo 

|e^t{«-eo) / rf^^ -Q F,((e.),ey(T;.)) 

(2.13) 

hence the result. □ 

Skeleton integrals are the fundamental objects from which regularized 
rough paths will be constructed in the next subsections. 
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2.2 Coproduct structure and increment-boundeiry decompo- 
sition 

Consider for an example the trunk tree T^'^" (see subsection 2.4 for an ex- 
planation of the notation) with vertices n^n — and labels 
i : {1, . . . , n} {1, . . . ,d}, and the associated iterated integral (assuming 
r = (r(l), . . . , r{d)) is a smooth path) 

ft rXn-l 

[lins = [/Tw„(r)]t. = / dr,M^)) ... dr,Mn)). (2.i4) 

J s J s 

Cutting T^"^" at some vertex v G {2, . . . , n} produces two trees, L^jT^^^" 
{left or rather bottom part of T^^") and i?^T^^" {right or top part), with 
respective vertex subsets {1, . . . ,v — 1} and {v, . . . ,n}. One should actually 
see the couple {L^T^'^" , R^T^'^") as L,,.T^'^" (g, R^T^'^" sitting in the tensor 
product algebra T(SiT. Then multiplicative property (ii) in the Introduction 
reads 

[SIjId„{T)]tus = ^ [-^L„TMn(r)]to[-fil^TWn(r)]„s. (2-15) 

„6y(TMn)\{i} 

On the other hand, one may rewrite [IjUn {T)]ts as the sum of the incre- 
ment term 

[6G]ts = /* dTMi)) r rfr.,(^(2)) . . . r-- dr^Mn)) 

- r dv^dm) r dr.Am ■ ■ ■ r-' dv^Mn)) 

(2.16) 

and of the boundary term 

dVM^))... dr,„^(£(ni)). 

/S rXn-i + l rXn-l 

dr^^^^Mni + 1)) J dV^^^^Mni + 2)) . . . y dr,Mn)). 

(2.17) 

The above decomposition is fairly obvious for n = 2 (see Introduction) 
and obtained by easy induction for general n. Thus (using tree notation this 
time) 

[Ilid.n{T)]ts = [5SkIjid„]ts — ^ [-fL„TWn(r)]ts ■ [Sk/^TWn(^)]s• 

^;eV(TWn)\{l} 

(2.18) 
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The above considerations extend to arbitrary trees (or also forests) as 
follows. 

Definition 2.8 (admissible cuts) 1. Let T be a tree, with set of ver- 
tices V(T) and root denoted by 0. If v = (vi, . . . ,vj), J > 1 is 
any totally disconnected subset of V{T) \ {0}, i.e. v-i -/» Vj for all 
i,j = 1, . . . ,J, then we shall say that v is an admissible cut ofT, and 
write V \= V{T). We let R^T be the sub-forest (or sub-tree if J = 1) 
obtained by keeping only the vertices above v, i.e. y(ii^T) = vyj{w G 
y(T) : 3j = 1, . . . , J, If -» Vj}, and L^T be the sub-tree obtained by 
keeping all other vertices. 

2. Let T = Ti . . . T; be a forest, together with its decomposition into trees. 
Then an admissible cut of T is a disjoint union ViL) . . .D Vi, Vi C Tj, 
where Vi is either 0, {Oj} (root ofTi) or an admissible cut o/Tj. By 
definition, we let L^T = L^^Ti . . . L^iTi, R^T = Rv^Ti . . . Rvi'^i (if 
Vi = 0, resp. {OJ, then (L^Ti,i2„Ti) := (Ti,0), resp. {%,Ti)). 

We exclude by convention the two trivial cuts U . . . U and {Oi} U 

...u{oj. 

See Fig. [2] and [3l Defining the co-product operation A : T ^ T ®T , 
T — > e®T + T®e + ^.y^y(Y>) L^T R^T (where e stands for the empty tree, 
which is the unit of the algebra) yields a coalgebra structure on T which 
makes it (once the antipode - which we do not need here ~ is defined) a Hopf 
algebra (see articles by A. Connes and D. Kreimer [3 [3 [8]). The convention 
is usuall to write v = c (cut), L^T = i?^(T) (root part), R^T = P'^(T) and 
A(T) = e®T + T®e + Y,c P'^i^) ® ^''i'^) (^ot^ the inversion of the order 
of the factors in the tensor product). 



w' 




Figure 2: Admissible cut. 
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Figure 3: Non-admissible cut. 



Eq. ()2.15j) extends to the general formula (called: tree multiplicative 
property), which one can find in [21] or |17j . 

[5IjiT)]tus= Yl [hMntuilRMnus, (2.19) 

v\=V(T) 

satisfied by any regular path F for any tree T. 
Letting formally s = itioo in eq. (j2.19p yields 



Mr)]tu = [sskhitu - i^LM^)]tu . [Sk/R„T(F)]n. (2.20) 

feV(T)\{o} 

which generalizes eq. (j2.18p . Conversely, eq. (j2.20p implies the tree 
multiplicative property eq. (j2.19p . as shown in Lemma 12.101 below. 

2.3 Regularization procedure 

Definition 2.9 (regularization procedure for skeleton integrals) Let 

T = {vi < . . . < v^f^} be a tree, T C T a subtree, fj, a compactly supported 

Borel measure on such that supp/i C ■ ■ ■ ,£,\y(f-)\) \ |^i| < ... < 

l^|y(T)|}' '^'^'^ Dreg C a Borel subset. 

The (formal) D^eg -regularized skeleton integral TlSklj is the linear map- 
ping (see eq. i2.9\) ) 

^ [7^SklT(M)]^ = (2^)-l^mi/2(;i,lD..,(6 • [sklT {{x.).evm - e^^--m-^5^) 

(2.21) 

where fl is the partial Fourier transform of fi with respect to ixv)veV(T)- 

By assumption we shall only allow Dreg = M z/ T is a tree reduced to one 
vertex. 
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Lemma 2.10 (regularization) Let T = Ti . . . be a forest, together with 
its tree decomposition. Define by induction on \ V{T)\ the regularized inte- 
gration operator [TZIj]^^ by 

il\[6nSkh^l^- J2 [^^^.T,],J^SkIfi,T,] i (2.22) 
Then [TZIj]^^ satisfies the following tree multiplicative property: 

[^'T^h],^s= E . [7^/i^„T]„,. (2.23) 

v\=V(J) 

By analogy with eq. [Km\2Al\WAm . [^7^SkITJ^^, resp. [7^/T^. (5)]t^ := 



~ 'Ylv\=v{T) \^^L^Tj\g \^^^'^Rvfj^ g fniy be called t/ie increment, resp. bound- 
ary operators associated to the tree Tj. 

Remark. By Definition 12.91 the condition [T^/flts = [ItIis holds for a 
tree reduced to one vertex. This implies in the end that one has constructed 
a rough path over the original path F. 

Proof. If the multiplicative property (|2.23|) holds for trees, then it 
holds automatically for forests since ['7^/Ti...Tjts is the product nj=i[^-^Tj]ts- 
Hence we may assume that T is a tree, say, with n vertices. Suppose (by 
induction) that the above multiplicative property (j2.23p holds for all trees 
with < n — 1 vertices. Then 

= E E (- [^^i-oi«m]j„ ['^■^^^^oL„(T)]„J7^SkI/^^T]s 

(2.24) 

Let X = V ]1 w := V \J w \ {i ^ v \J w \ 3j ^ v \J w \ i-» j}. Then 
one easily proves that o L^iT) = L^iT), R^i^) = Rv ° ^£c(T) and 
Rw ° Lvi^) = L^o Rx(T). Hence 

[dTZhjtus = E [^^^="T]tJ- E ['^^L^{R.TMnskiR^^R^j^]s + mskiR^j], 

x\=V{T) \ v\=V(Ra,T) 

= E [T^lL^T]tu[niR^j]us. (2.25) 

x\=ViT) 

□ 
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2.4 Permutation graphs 

Consider now a permutation a G S^,. Applying Fubini's theorem yields 

J s ^ s ^ s 

rti rt'i rtn 

= / dr,^(,,(%(l))) / dr,^(,,(%(2)))... / dT,^^Ji{a{n))), 

J Sl J S2 J Sn 

(2.26) 

with si = s,ti=t and Sj E {s} U i < j}, tj G {t} U i < j} (j > 

2). Now decompose dr,^^^, (£(a(i))) into -/;^) dr,^,^,, (^(a(j))) if 

/ s,t, + t, and /^*.dr.^(^,,(^(a(j))) into (/J - /;^) dr.^^^,^ (^(a(j))) if 
Sj ^ s. Then /^(r) has been rewritten as a sum of terms of the form 

± r (ir,,(^(a(l))) r dr^M^{2))) . . . r dV^Mc^in))), (2.27) 

t/ S i/ S " s 

where ti = t and Tj S {t} U {xj, i < j}, j = 2, . . . ,n. Note the renaming of 
variables and vertices from eq. (j2.26p to eq. (|2.27p . Encoding each of these 
expressions by the forest T with set of vertices ^(T) = label 
function £oa, roots {j = 1, . . . ,n \ tj = t}, and oriented edges \ j = 

2, . . . , n, Tj = Xj-}, yields 

/^(r) = /T^(r) (2.28) 

for some T'^ G T called permutation graph associated to a. 

Summarizing: 

Lemma 2.11 (permutation graphs) To every permutation a ^ Tin is 
associated a permutation graph 

r = Y,gia,j)T'^ gT, (2.29) 
i=i 

9{(^ij) = if) each forest TJ being provided by construction with a total 
ordering compatible with its tree structure, image of the ordering {vi < 
. . . < Vn} of the trunk tree T^'^" by the permutation a. The label function of 
T°" is i o a, where i is the original label function of T^'^" . 
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Example 2.12 Let a = (^^ ^ ^ j . Then 

f dV^M^)) r dT^M'^)) r dT^M^)) = 

J s J s J s 

j-t j-X2 I-X2 

- / dr,,(^(2)) / dT^M^)) / dVM^)) 

+ / dr,2(^(2)) / dv^M^)) ■ / dr,i(^(i))- (2.30) 

OS iJ s ^ s 

Hence T'^ = — + T2 is the sum of a tree and of a forest with two compo- 
nents (see Fig. 



3 2 



1 1 



3 



2 3 



Figure 4: ExampleEH From left to right: T^,TJ; L{i}T^ «> i?{i}T^; L{i.2}T^ 

^{1.2}^ 



2.5 Fourier normal ordering algorithm 

Let r = (r(l), . . . , r((i)) be a compactly supported, smooth path, and 
r"(zi, . . . , i„) some iterated integral of T. To regularize V^{ii, . . . , we 
shall apply the following algorithm (a priori formal, since skeleton integrals 
may be infra-red divergent) : 

1. (Fourier projections) Split the measure /U = dT{ii) ... dT{in) into 
T.a^ii^^-H'^D^miO). where D'^ = G | 1^,(1) | < 

■ ■ ■ ^ |^CT(n)|}) s-iid /U is the Fourier transform of fi. We shall write 



:= {Id" -fj) ° <y = (l£,id„.(/i o fj)) ; 



(2.311 



2. Rewrite li (.F-I(Id-A)) , where ^(j) = i„ as /t^Ca^'^) := E/=i 5(fT, jVtj (a^"), 
where T°" is the permutation graph defined in subsection 2.4; 
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3. Replace If a (/i'^) with some regularized integral as in Definition 12.91 and 
Lemma 12.101 

nij^ifi") ■.= ^gia,j)nijjifin; (2.32) 

4. Sum the terms corresponding to all possible permutations, yielding 
ultimately 

7er"(n,...,i„) = 7^/T.(^'^). (2.33) 

Explicit formulas for T = B"^ may be found in the following section. 
Theorem 2.1 

IZV satisfies the multiplicative (ii) and geometric (Hi) properties defined 
in the Introduction. 

The proof given in [36] shows actually that any choice of linear maps 
[T^Sk/xjf : pL — > [7^Sk/x(/u)]t such that 

(i) [7^Sk/Tl.T2(/^l ® l^2)]t = [■7^Sk/Tl(/Ul)]^[7^Sk/T2(/U2)]t and 

(ii) [7^Sk/T(/)]^ = [Sk/T(/)]t = J* f{u) du if T is the trivial tree with 
one vertex, 

yields a regularized rough path over F if F is smooth. Hence our 'cut' 
Fourier domain construction is arbitrary if convenient. As already said in 
the Introduction, it seems natural to look for some more restrictive rules 
for the regularization; iterated renormalization schemes (such as BPHZ or 
dimensional regularization) are obvious candidates (work in progress). The 
question is: is such or such regularization scheme better in any sense ? Con- 
trary to the case of quantum field theory where all renormalization schemes 
may be implemented by local counterterms, which amount to a change of the 
value of the (finite number of) parameters in the functional integral (which 
are experimentally measurable), and give ultimately after resumming the 
perturbation series one and only one theory, we do not know of any proh- 
abilistically motivated reason to choose a particular regularization scheme 
here. 

3 Rough path construction for fBm: case of dis- 
tinct indices 

The strategy is now to choose an appropriate regularization procedure, so 
that regularized skeleton integrals of are finite and satisfy the uniform 
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Holder and convergence rate estimates given in Theorem 10. 1[ 

3.1 Analytic approximation of fBm 

Recall B may be defined via the harmonizable representation [32] 

Bt = Co. [ m"^"'^—^WidO (3.1) 

where {W^,(, > 0) is a complex Brownian motion extendeded to M by setting 
W.^ = -W^ > 0), and c« = 1^33=^. 

We shall use the following approximation of i? by a family of centered 
Gaussian processes {B^i,ri > 0) living in the first chaos of B. 

Definition 3.1 (approximation B"^) Let, for r] > 0, 

B'l = c^f e-^l«l|e|^-°^^^ W{di). (3.2) 

The process B"^ is easily seen to have a.s. smooth paths. The infinites- 
imal covariance E(5^)^(iJ^)J may be computed explicitly using the Fourier 
transform [llj 

.F<-(0 = -|= / K'^-{x)e-'^^dx = — ^^^e-^^l^licr-^.^ 

V 27r Jr ' 2cos7rar(-2a) 

(3.3) 

where {s — t) := 2ros^"'' (~i('^ ~ + 2r/)^°~^. By taking the real part 
of these expressions, one finds that B^ has the same law as the analytic 
approximation of i? defined in [34], namely, B^ = rj_|_i^ + rt_i^ = 2Re r^+i^, 
where F is the analytic fractional Brownian motion (see also [33j). 

3.2 Choice of the regularization procedure 

Let a E S„ be a permutation. Recall (see Lemma [2. lip that the permutation 
graph T'^ may be written as a finite sum 'Yl,'j=i9{'^i ^ where each TJ is 
a forest which is automatically provided with a total ordering. In the two 
following subsections, we shall consider regularized tree or skeleton integrals, 
IZIf or T^Sk/x, for a forest T which is one of the TJ. 

Definition 3.2 Fix Creg £ (0)1)- Let, for T with set of vertices ViT) = 

{Vi< ... < Vj}, 

:={(e.i,...,e.,)GM^ I |e.J<...<|e.,|}, (3.4) 
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^Jeg ■■= • • • ,^.,) G I Vz; G V(T), \Cv+Y, U > C.eg max{|e^|; w ^ 

(3.5) 

andlZIj, resp. TZSklj be the corresponding Mj^g-regularized iterated, resp. 
skeleton integrals as in subsection 2.3. 

Condition (j3.5p ensures that the denominators in the skeleton integrals 
are not too small (see Lemma l2.7p . 

The following Lemma (close to arguments used in the study of random 
Fourier series [20j) is fundamental for the estimates of the following subsec- 
tions. 

Lemma 3.3 (i) Let F{u) = j^dW^a{i)e^< , where |a(OP < C|Cr^~^^ 
for some < /3 < !.■ then, for every ui,U2 G M, 

E|F(tii) - F(ti2)|2 < C'\ui - U2\^^. (3.6) 

(ii) Let F{r]) = ^ dW^a(Oe-^l^l (r? > 0), where |a(OP < C\^\-^-^^ for 
some < l3 < 1: then, for every 7/1,772 G M+, 

E|F(r/i)-F(7?2)|2<C7'|7?i-7?2p^. (3.7) 

Proof. Bound |e^"i« - e'"25| by |7ii - U2||^| for |^| < and by 2 

otherwise, and similarly for |e~''^l^l — e~''2l^l|. Note the variance integral is 
infra-red convergent near = 0. □ 
Remark: Unless |a(C)P is Lj^^ near ^ = 0, only the increments F{ui) — 
F{u2), F{rji) — F{rj2) are well-defined. 

3.3 Estimates for the increment term 

In this paragraph, as in the next one, we consider regularized tree integrals 
associated to 7^B"'''(ii, . . . ,i„,) where ii ^ . . . ^ in are distinct indices, so 
that B(ii), . . . , B{in) are independent. 

Lemma 3.4 (Holder estimate and rate of convergence) Let T = TJ 

for some j, and a < 1/\V {T)\. 

1. The skeleton term 

[G^''^(il, . . .,in)]u = [7^SkIT {{dB'^iii) ... dB^{in)Y)]^ (3.8) 
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(see eq. 112.31]) ) writes 



satisfies the uniform Holder estimate: 



E\[6G:i;''{h, . . .,in)]tsf < C\t - (3.10) 

3. (rate of convergence) : there exists a constant C > such that, for 
every r]i,rj2 > and s,t G M, 

E I [SGl^^'^iiu- ■ ■ , in)]ts - . . . , in)]ts\^ < C\m - ml''"- 

(3.11) 

Proof. 

1. Follows from Lemma 12.71 and the definitions of B"^ and of regularized 
integrals in the previous subsections 2.3 and 3.1. 

2. (Holder estimate) 

One may just as well (by multiplying the integral estimates on each 
tree component) assume T is a tree, i.e. T is connected. 

Let V{T) = {vi < ... < v\v{j)\}, so that \Cv^ \ < ... < |C«|v(t)|I- Since 
every vertex v E V{T) \{vi} connects to the root vi, one has 

I^C^)! • l?t'|y(T)|l — + • • • + ^'i'|XA(T)| I ^ G»'egl'?-y|v(T)| I ' (3.12) 

so that ^ := J2veV{J)^'" comparable to £,v^y(j-j\, i-e. belongs to 
[C~'^Cv\Y(j)\ 5 G^''|v(T)|] some large enough positive constant. Write 

[Gr(n,---,^n)]n = /Ke^"Me)de 

Vertices at which 2 or more branches join are called nodes, and vertices 
to which no vertex is connected are called leaves (see Fig. [5]). 

The set Br{vi — » ^2) of vertices from a leaf or a node vi to a node 
V2 (or to the root) is called a branch if it does not contain any other 
node. By convention, Br{vi V2) includes vi and excludes V2- 
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Figure 5: 3,4,6 are leaves; 1, 2 and 5 are nodes, 2 and 5 are uppermost; branches 
are e.g. Br{2 1) or Br{6 ~» 1). 

Consider an uppermost node n, i.e. a node to which no other node is 
connected, together with the set of leaves {wi < . . . < wj} above n. 



Let pj = \V{Br{wj n))\. Note that — 7 



1 \ 2 



< 



Wj to n 



Now we proceed to estimate Var a(^). On the branch number j from 



\iv\<\iwA,v&Br(wj^mi)\{wj} 



n 



v(iBr(wj — »n) 



< 



-l—2apj 



(3.13) 



and (summing over . . . , Cwj^i and over ^n) 



lf-j_il<IC-jl 

fjf If \-l-2api 



-1— 2apj_ 



l«»ll<l«»2l 
-(l+2apj)+[2-2a(l+pi+...+p.j-i)]-2 



l-2a 



ieni<ic»ii 



Wj 



-l-2aW{n) 



(3 



where W{n) = pi + ...+ pj + 1 = \{v : v ^ n}\ + 1 is the weight of n. 

One may then consider the reduced tree T„ obtained by shrinking 
all vertices above n (including n) to one vertex with weight W{n) 
and perform the same operations on T„. Repeat this inductively 
until T is shrunk to one point. In the end, one gets Var a(^) < 
l^^lvmiT^"^"'^^^^' ^ ^^^"1^(^)1. Now apply Lemma EJ (i) . 



3. (rate of convergence) 
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Let X^""' := [G^i''^(ii,...,i„)]„ - Expanding 



\V(T)\ 



^ g-'?2(ie«il+-+|^_J)(-g-ml^l _ g-r?2|^|^g-m(l^+il+-+l^.v-{T)l) 
i=i 



gives X:?^'"^ as a sum, X^I^'^^ = E.ey{T) X^'^'^'iv), where 
/de.^4ef)(e-^^'^"' _e-^2ie.|) is obtained from [G^ {h, ■ ■ ■ ,in)]u by 
replacing e"^!^"! with e"^!!^"! — e"''^!?!.!^ and e~''l^™l, w ^ v either by 
g-»?ilCiu| Qp )-,y g-f)2|Cm|_ want to estimate Var bu{(,v) uniformly in u. 

Fix the value of .^^ in the computations in the above proof for the 
Holder estimate. Let wj be the maximal leaf above v, and n ^ be 
the node just above f if u is not a node, n = v otherwise. Summing 
over all nodes above n and taking the variance leads to an expression 
bounded by \iwj\~^~'^"''^^'^\ where W{n) = \{w : w ^ n}| + 1 is 
as before the weight of n. Consider now the corresponding shrunk 
tree T„. Let Tn{v) be the trunk tree defined by Tn{v) = {if G T„ : 
w ^ V oi V ^ yj {v}; similarly, let T{v) be the tree defined by 
T{v) = {w ^ T : w ^ V OY V ^ {u}, so that Tn{v) is the 

corresponding shrunk tree. Sum over all vertices w E T„(f) The 
variance of the coefficient of e^''^!^"' is 



< 



i€ni>ie.i 



dS,n\S,n\ 



-l~2aW{n) 



-l-2a 



\in,\<\iu\,weJ„^{v)\{n,v} 



< 



w<^Tn(v)\{n,v} 



(3.15) 



if u 7^ n, and 



IC»l<l5n|,«'6T„(^)\{n} ^gTr„(i,)\{n} 

(3.16) 



if V 



n. 



Removing the vertices belonging to T{v) from T leads to a forest which 
gives a finite contribution to the variance. Hence (by Lemma 13.31 (ii)) 
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□ 

The notion of weight W{v) of a vertex v introduced in this proof wih be 
used again in subsections 3.4 and 4.1. 

3.4 Estimates for boundary terms 

Let T = TJ for some a G and ii / . . . 7^ z„ as in the previous subsection. 
By multiplying the estimates on each tree component, one may just as well 
assume T is a tree, i.e. is connected. 

We shall now prove estimates for the boundary term TZIj ({dB^{ii) ® . . . dB'^{in)Y) (d) 
associated to T (see Lemma l2.10p . 

Lemma 3.5 Let T = TJ for some j (so that n = \V{T)\). 

1. (Holder estimate) The regularized boundary termTZIj {{dB^{ii) ... dB^{in)Y) (d) 
satisfies: 

E |[7^/T iidB^iii) ® . . . C5 dB^iinW) {d)],^ p < sp^l^WI (3.17) 
for a certain constant C. 

2. ( rate of convergence ) There exists a positive constant C such that, for 
every r/i, 772 > 0, 

E|[7^/T ((dS^Hn) ® . . . ® dB^'iinW) {d)]ts- 

-[Uh {{dB^Hk) ® . . . ® dB^'{in)r) mts\^ < c\m - ml''"- 

(3.18) 



Proof. 

1. Apply repeatedly LemmaEIOlto T: in the end, [TZh {{dB'^{ii) ® . . . dB'^{in)T) {d)]ts 
appears as a sum of 'skeleton-type' terms of the form (see Figure E]) 

Ats := [6nSklLT]ts ■ [7^SkIR^^„i^^_^„...<,J^^^(T)]s... [7^SkI^^^„i^^(T)]s[7^SkIK^^^]^ 

{{dB'^{ii)^...(^dB'^{in)T), 

(3.19) 

where vi = {vi^i < ... < w^jj |= T, V2 |= L^^T, vi = {vi^i < 
... < viji^) 1= L^,;_^ o . . . o L^,j(T) and LT := o . . . o L„^(T). In 
eq. ()3.19p the forest T has been split into a number of sub-forests. 
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LTD (Uj^^Tjj ; we call this splitting the splitting associated to A^g for 
further reference. 

First step. 

Let [B^' [^]]u n/ii be the contribution to TZSklR^ j of 

all Fourier components such that ^ = {Cvi^i, ■ ■ ■ jCvu^^), \^vi,i\ < • • • < 
l^^j J J is fixed. For definiteness (see Definition 13. 2p . 



Ji 



(3.20) 



Then 



VarKMC]]. 



hence 



n 

1 



-1-2q 



n i^-r'"'" 

(3.21) 



(3.22) 



Second step. 



More generally, let ^^^^^ [^] H/Li '^^5^; ■ j)) be the contribution 
to 

[nSkln^^,L.^_^,,„,L.^^j)l . . . [7^SkIH,^oL„,(T)].[^SkI^„^,], iidB^iii) g 

(3.23) 

of all Fourier components such that ^ = (^^,; ^ , . . . , ^t,; ) is fixed. Then 



Var(i?J--'-q^]) < n le.r'l^^^"^)!"-^ 



(3.24) 



(proof by induction on /). 
Third step. 



dB^inW) 
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Let V{LT) = {wi < ... < Wmax}- By definition, Ats = as(H)(e''='* ■ 
e^^^)dE, with 

as{E)= [d^ [... [ n dW^J^y^)) 



n.ey(LT)(-iCa)e-''l«-l|e»l^-" 



(3.25) 

where Fourier components in satisfy in particular the following 
conditions: 

• + J2w'^w,w'&V{LJ)^w'\ > Cregmax{|^^/| : w' w,w' e 
V(LT)}; in particular, ( , l^"'^'" ^ ] < ICwt^'^"] 

• Z]«,ey{LT) = 

• for every u; G F(LT), < and < for every 
V £ R{w) := {v = vi^i, . . . ,vij^ \ V ^ w} (note that R{w) may 
be empty). See Fig. [H 

2 (• .0 "2,1 



Figure 6: Here V{LT) = {0,1,2,4}, R{0) = i?(4) = 0, R(l) = {w2,i}, i?(2) = 

{V2.2}- 

Note that |H| < |C«;maa;l ~ 1^1 since every vertex in V{LT) connects to 
the root (see first lines of the proof of Lemma 13.41 (2)). 

If w G LT, split R{w) into /?(?/;)> U R{w)<:, where R{w)^ := {v £ 
R{w) I V ^ Wmax}- Summing over indices corresponding to vertices in 
i?T> := {v = vi^i, . . .,vij^ I V > Wmax} = Uyj(zLTRiw)> , One gets (see 
again proof of Lemma 13.41 (2)) 

-Q f ^ ^^^|^^|-2|y(R„T)|a-l < |^|-2aE„e«T>inR.T)|_ ^3^26) 
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Let w £ LT\{'Wmax} such that i?(w)< 7^ (note that R{wmax)< = O)- 
Let R{w)<: = {vi-^ < ... < Vi-}. Then (integrating over {(,v),v £ 



J\iv,^\>\u J\iv,^\>\iv,^\ J\iv,.\>\iv,^_^\ 

(3.27) 

In other words, each vertex w G LT 'behaves' as if it had a weight 
1 + X]i;ei?(w))< l^(-^t>1')l- Hence (by the same method as in the proof of 

Lemma[3a](2)) Var(a,(0) < |Hr^-^'^(l^(")l+^--< I^^^^^^D . iSf^"^— > 1^^^ 
|^|_i_2a|y(T)|^ Now apply Lemma [23] (i) . 

2. Similar to the proof of Lemma 13.41 (3). Details are left to the reader. 

□ 



4 End of proof and final remarks 

4.1 Estimates: case of coinciding indices 

Our previous estimates for E|7^B"j''(zi, . . . ,in)p (Holder estimate) and 
E|7^B"j''^ (zi, . . . ,i„) — 7^B"j^^(zi, . . . ,in)\^ (rate of convergence) with ii ^ 
. . . ^ in rest on the independence of the Brownian motions W{ii), . . . , VF(i„). 
We claim that the same estimates also hold true for E|7?.B"''''(zi, . . . ,in)P 
and KlTZB^f^ {ii, . . . , in)— TZB^f^ {ii, . . . , in)\'^ if some of the indices (ii, . . . , in) 
coincide, with the same definition of the regularization procedure TZ. The 
key Lemma for the proof is 

Lemma 4.1 (Wick's lemma) (see f2B], ^5.1.2 and 9.3.4) 

Let {Xi, . . . , Xn) be a centered Gaussian vector. Denote by Xi-^o. . -oXj^. 
(1 < ii , . . . , ifc < n) or : X-i^ . . . Xi^ : the Wick product of Xi^ , . . . , Xi^ ( also 
called: normal ordering of the product Xj^ . . . Xi^_), i.e. the projection of the 
product Xi^ . . . Xi^ onto the k-th chaos of the Gaussian space generated by 
Xi, . . . , Xn. Then: 
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1. 



Xi . . . X„ = o . . . o X„ + ^ ¥.[Xi^Xi.,]Xi o ...oXi^o ...oXi^o ...oXn 



(il.«2) 

^[Xi^Xi^] . . . ¥.[Xi^^^^Xi^^^^] 



(n,«2),.--:(«2fc+l,»2fc+2) 

Xi o . . . o o . . . o o . . . o Xj2^_^i o . . . o o . . . o X„ 

+ ..., (4.1) 

where the sum ranges over all partial pairings of indices {11,12)-, ■■■ ■, {i2k+iii2k+2) 
(l<A;<LfJ-l). 

2. For every set of indices ii, . . . ,ij,i'i, . . . 



E 



(X,,o...oX,^,)(X,,o...oX,,)j = Il^i^^mXi' ]. (4.2) 



In our case (considering Tl'B^^^{ii, . . . , i„)) we get a decomposition of the 
product dW^^{ii) . . . dW^^{in) into (iWg^(ii) o . . . o dW^^{in), plus the sum 
over all possible non-trivial pair contractions, schematically {W^,{ij)W^^{ij')) = 

Consider first the normal ordering of TZB^^^{ii, . . . , i„). As in the proof 
of Lemma 5.10 in [33], let be the 'index-fixing' subgroup of S„ such that 
: 0"' E <^=^> Vj = 1, . . . ,n, = ij- Then (by Wick's lemma and the 

Cauchy-Schwarz inequality) : 

Var :7^BS^(il,...,i„) := E |: 7^B^^;''(n, . . . , i„) :|' 

= E [: 7^BS''(1, . . . , n) : : 7^BS''(a'(l), . . . , a'(n)) :] 

< . E|7^B"•''(l,...,n)|^ (4.3) 

hence the Holder and rate estimates of section 3 also hold for 
:7^B-•'?(il,...,i„) :. 

One must now prove that the estimates of section 3 hold true for all pos- 
sible contractions of 7^B"'^(zi, . . . ,in)- Fixing some non-trivial contraction 
(jii J2)i • • • 5 (i2Z-ii J2«)' ^ ^ li results in an expression X™"*'" belonging to 
the chaos of order n — 21. By necessity, ij^ = ij^, . . . ,ij2;-i ~ ^21^ 
may well be that there are other index coincidences. The same reasoning 
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as in the case of : TZ'B^f'{ii, . . . ,in) '■ (see eq. ()4.3p ) shows that one may 
actually assume im / im' if m / m' and {m,m'} / {ji,j2}, ■ ■ ■ ,{j2i~ij2i}- 
Now (as we shall presently prove) the tree integrals related to the contracted 
iterated integral Xj""**" may be estimated by considering the tree integrals 
related to Xj,, := 7^B"g~^''''(ii, . . . . . . , ij^j , . . . ,in) (which has same law 

as 7^B"g~^'''^(l, . . . , n — 2/)) and (following the idea introduced in the course 



of the proof of Lemma 13. 4p increasing by one the weight W of some other 
(possibly coinciding) indices j'l, ■ ■ ■ ,j'2i / ji, • • • ,j2i ~ or, in other words, 'in- 
serting' a factor . . . |Cj^,J~^" i^i the variance integrals -. This amounts 
in the end to increasing the Holder regularity {n—2l)a^ of X^^ by 2/a, which 
gives the expected regularity. 

Fix some permutation a £ Tin, and consider the integral over the Fourier 
domain |^o-(i)| < ... < |Co-(n)l as in section 2. Change as before the order 
of integration and the names of the indices so that dW^^^^j (ij) — > dW^^ (V{j)); 
for convenience, we shall still index the pairing indices as (ji , j2 ) , . . . , {j2i- 1 , j2i ) ■ 
We may assume that \j2k-1 — j2k\ = 1, k = 1,...,/ (otherwise = 
I'^i2fe-il = foi^ hk-i <m < j2k or j2k <m < j2k-i, which corresponds 

to a Fourier subdomain of zero Lebesgue measure). In the sequel, we fix 
a and (j,/) = (j2A:-i, J2fc) for some k. 

Let T = Ti...Tl be a forest appearing in the decomposition of the 
permutation graph T*^ as in subsection 2.4. Applying repeatedly Lemma 
12.101 to T leads to a sum of terms obtained from the contraction of Ats = 
Ats{l)...Ats{L),withAtsik) = [SnSkl^fJts UjlT^SklrJs (^((S),^y(^f^/B^ii,: 

where LTi~, T'^ . . . , T'^ . . . are all subtrees appearing in the splitting as- 
sociated to Ats{k) (see proof of Lemma [33]) . 

Let T be one of the above trees, either LT/j or T'^ j. Reconsider the proof 
of the Holder estimate or rate of convergence in 
The integrals 



Skl((x,),ev(T) ^e^^-evm-"?'' 



emma 13^ or Lemma 
appearing in the defini- 
tion of the regularized skeleton integrals write i"!^^'''')! -n — - — /T^^ t ^^ 

Llvev(j)\t,v+2^^u^^ l;w) 

(see Lemma 12. 7p . After the contractions, one must sum over Fourier in- 
dices {Cv)v(^v{T) such that (Ci.)t,eV(T) G ^Jeg and Ci2m-i = "62™ if both 

J2m-1, j2m S ^(T). 

Let T be the contracted tree obtained by 'skipping' {ji, . . . , ^2/} H V(T) 
while going down the tree T (see Fig. [71 [H [9]). 

The denominator ICv+XltoeT w^v is larger (up to a constant) than the 
denominator + "l^w&f w^v^'>"\ obtained by considering the same term in 
the contracted tree integral X^^ (namely, ~^Ylw€T w^v^^l is of the same 
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Figure 7: Case (i-a). T and f . 




Figure 8: Case (i-b). T and f . 



order as max{|^^|; tt; £ T,w ^ v} > max{|^^|;ii; G T, tt; -» v}). Hence 
E(A^°"'*'')^ may be bounded in the same way as ^A^^ in the proof of Lemma 
13.41 or Lemma [331 except that each term in the sum over {S,v,v G V{T),v / 
ii)---)j2i) comes with an extra multiphcative pre-factor 5 = S{{£,y),v £ 
V{T),v / jl, . . . ,j2i) - due to the sum over {Cjm)m=i,...,2i ^ which may be 
seen as an 'insertion'. 

Let us estimate this prefactor. We shall assume for the sake of clarity 
that there is a single contraction (ii,j2) = U^j') (otherwise the prefactor 
should be evaluated by contracting each tree in several stages, 'skipping' 
successively (ji, ^2), ■ ■ ■ , {j2i-i,j2i) by pairs). As already mentioned before, 
\j — j'l = 1 so that j and / must be successive vertices if they belong to 
the same branch of the same tree T. Note that, if j and / are on the same 
tree, the Fourier index H := YlveV{T) (used in the Fourier decomposition 
of Lemma 13.41 or in the third step of Lemma 13. 5p is left unchanged since 

e, + e,' = o. 

Case (i): (j, j') belong to unconnected branches of the same tree T. This 
case splits into three different subcases: 

(i-a) neither j nor j' is a leaf. Let w, resp. w' be the leaf above j, resp. j' 
of maximal index and assume (without loss of generality) that ^ 
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Figure 9: Case (i-c). T and f . 



le^'l- Then 



l—2a \ 



(4.4) 

which has the effect of increasing the weight W{w) by 2. 

(i-b) j is a leaf, j' is not. Let w' be the leaf of maximal index above /. 
Then 



(i-c) both j and j' are leaves. Let v, resp. f' be the vertex below j, resp. 
/, i.e. j V, j' — > v' . Then 

5<f/" dC,\^,\-'-A <m-'^ (4.6) 

\J|5j|>max(|C„|,|OI) / 

which has the effect of increasing W{v) by 2. 

Case (ii): are successive vertices on the same branch of the same 

tree T. Assume (without loss of generality) that j — > j' . Then S" = if j is 
a leaf (since ^j/ + Ylw^j' = + Cj' = and such indices fail to meet the 
condition defining MjFg^), otherwise S < if w is the leaf of maximal 

index above j (by the same argument as in case (i-a)). 

Case (iii): belong to two different trees, T and T'. 

This case is a variant of case (i). Nothing changes compared to case (i) 
unless (as in the proof of Lemma 13.41 or in the 3rd step of Lemma 13. 5|) one 



(4 
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needs to compute the variance of the coefficient a(H) or as(H) of e™" for H 
fixed. Assume j belongs to the tree T = LT^ while j' is on one of the cut 
trees T'^_^, . . . , T'^^^-, . . . 

Assume first j is not a leaf, and let w be the leaf above j. Then the 
presence of the extra vertex j modifies the Fourier index H in the Fourier 
decomposition of A^°'^^'' (k) , = /r a(S)(e'^*-e'^')dS or Afr*''!^) = 

Jjg a5(H)(e'"* — e'"*)(iH, by a factor which is bounded and bounded away from 
0, hence S < l^^u;]"^" as in case (i-a). 

If j is a leaf as in case (i-b) - while w' is as before the leaf of maximal 
index over j' -, one has: < |H| < Hence the sum over contributes 
an extra multiplicative pre-factor S to the variance of the coefficient of a(H) 
or as(^) of order 

s<([ ^e,T^) <([ <l-r^ 



'|=|/2<|Cj|<2|H| \^j^w'\ J ~ \y|H|/2<|5,|<2|H| 

(4.7) 

which increases the Holder index by 2a (see Lemma l3.3p . 

The case when both j and j' belong to left parts LTf^., LT^/ is similar 
and left to the reader. □ 

This concludes at last the proof of Theorem 10. 1[ 

4.2 A remark: about the two-dimensional antisymmetric 
fBm 

Consider a one-dimensional analytic fractional Brownian motion F as in |33j . 

Definition 4.2 Let Zt = {Zt{l), Zt{2)) = (2Re Ft,2Im Ft), t e R. We 
call this new centered Gaussian process indexed by M the two-dimensional 
antisymmetric fBm. 

Its paths are a.s. Q~-H61der. The marginal processes Z{1), Z{2) are 
usual fractional Brownian motions. The covariance between Z{1) and Z{2) 
writes (see [33] ) 

Cov(Z,(l),Zt(2)) = -^f^[-sgn(s)|.|2"+sgn(t)|t|2"-sgn(t-s)|t-.|2°]. 

(4.8) 

Note that we never used any particular linear combination of the analytic/ anti- 
analytic components of B in the estimates of section 3 and 4. Hence these 
also hold for Z, which gives for free a rough path over Z satisfying Theorem 
10.11 of the Introduction. 
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